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ABSTRACT 

The Hubble Frontier Fields (HFF) are six clusters of galaxies, all showing indi¬ 
cations of recent mergers, which have recently been observed for lensed images. As 
such they are the natural laboratories to study the merging history of galaxy clusters. 
In this work, we explore the 2D power spectrum of the mass distribution Pyi(k) as a 
measure of substructure. We compare Pu{k) of these clusters (obtained using strong 
gravitational lensing) to that of ACDM simulated clusters of similar mass. To compute 
lensing Pm(M we produced free-form lensing mass reconstructions of HFF clusters, 
without any light traces mass (LTM) assumption. The inferred power at small scales 
tends to be larger if (i) the cluster is at lower redshift, and/or (ii) there are deeper 
observations and hence more lensed images. In contrast, lens reconstructions assum¬ 
ing LTM show higher power at small scales even with fewer lensed images; it appears 
the small scale power in the LTM reconstructions is dominated by light information, 
rather than the lensing data. The average lensing derived Pyi{k) shows lower power 
at small scales as compared to that of simulated clusters at redshift zero, both dark- 
matter only and hydrodynamical. The possible reasons are: (i) the available strong 
lensing data are limited in their effective spatial resolution on the mass distribution, 

(ii) HFF clusters have yet to build the small scale power they would have at z ~ 0, or 

(iii) simulations are somehow overestimating the small scale power. 

Key words: gravitational lensing: strong, galaxies: clusters: individual: Abell 
2744, Abell 370, Abell S1063, MACS J0416.1+2403, MACS J0717.5+3745, MACS 
J1149.5+2223 


1 INTRODUCTION 

Clusters of galaxies are the largest self-gravitating objects 
in the Universe. Their ultimate origins must he in some 
process of quantum fluctuations in an expanding uni verse 

fechrodingerlll939l : lHarrisorilll967l : lFedderke et al.l2015l ). but 

the earliest observable precursors of galaxy clusters are fluc¬ 
tuations in the cosmic microwave background. CMB fluctu¬ 
ations on the scale of individual clusters wo uld be at l ~ 10 4 , 
far in the diffusion-dampin g tail (iSilkll968i) and so far barely 
accessible observationally (iReichardt et al.l l2009). but their 
subsequent growth through linear gravitational instability is 
straightforward. 

* irshad@physik.uzh.ch 


The earl y _ observable proto-clusters (e.g, 

iToshikawa et al.1 120141 ) are, however, far beyond the 
baby fluctuations of the linear regime of gravitational insta¬ 
bility — in the well-known toy model of spherical collapse, 
the dynamical outer boundary of a cluster is at an overden¬ 
sity of ~ 200 times the critical density of the Universe. In 
this regime some analytical meth ods based on generalising 
spherical collapse are availabl e J Press &; Schechten Il974l : 

I Bond et al.l 1 199 it ISeliakl 120001 ) but, for the most part, 
theoretical study depends on numerical simulations. 

With the recent developments in computational re¬ 
sources, it is now possible to simulate dark matter and gas 
in cosmological volumes with good resolution. In such sim¬ 
ulations, it is possible to track the evolution and merg¬ 
ers of small systems to form large collapsed clusters of 
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galaxies. The agreement between simulations and observa¬ 
tions have been improving for various macroscopic proper- 
ties of galaxies, such as intergalactic g a s, bulge sizes etc . 


(Klvpin et al. 

2011: Angulo et al 

2012; Alimi et al. 2012j; 

Skillman et al. 

2014; Agertz et al. 

201l|: Voeelsberger et al.1 

20ll, 20l]l; [M 

arinacci et al. 2014 

: Schave et al. 20151). To 


study individual objects in mor e detail, hydrodynamical 
zoom-in simulations can be used (IFeldmann fe Mavedl2015l : 
iFiacconi et al.l 120151 ") in which individual objects are re¬ 
simulated at higher resolution using initial conditions de¬ 
rived from the larger-volume simulation. 

A more ambitious comparison between simulations and 
observations would be that in terms of the distribution of 
mass in the clusters. This is non-trivial for two reasons: 


• First, mass is not an observable; what we observe on the 
sky is light, mass can only be inferred. Tracing the mass typ¬ 
ically involves additional assumptions, such as that galaxies 
sit in the potential wells of the dark-matter, light traces mass 
with some scaling parameters, etc. 

• Second, since individual halos, galaxies or clusters of 
galaxies are the outcome of gravitational collapse and var¬ 
ious baryonic processes of a random held (initial density 
held), it is not possible to compare their mass distribution 
or clustering properties directly. The best one can hope for 
is to compare them statistically. 


The hrst problem could be solved with the help of grav¬ 
itational lensing, which is sensitive only to the total mass. 
Well-known examples of light not tracing mass, rev ealed 
by lensing, are the clusters ACO 2744 ([Me rten et a ll 1201 il l 
and recently ACO 3827 from Massev et al.l (120151 ') (see also 
I Williams fe Sahal 1201 ll : [Mohammed et al. 20141) . But simu¬ 
lations of cluster formation cannot be htted so as to repro¬ 
duce detailed properties of individual clusters. 

The second problem could be handled by identifying ro¬ 
bust statistical properties of the cl usters. The simp l est qua n- 
tity is the radial density prohles. iNewman et al) d2013bl fa|j 
studied the average density prohles of the lensing clusters 
and compared them to the simulated clusters. However, this 
quantity can only be measured if the cluster is virialised, 
dehnitely not for an ongoing merger. High-mass clusters in 
simulations do not usually virialise until z ~ 0.3; before that 
they often show elongation, multiple cores and many sub¬ 
structures indicating a recent merger, as do observed lensing 
clusters. 

Another very popular statistic is the mass function; 
countin g the number of subhalos as a function of their 
masses (Nataraian fe Springell 120041 ; iNataraian et al.l 120071 : 
lAtek et al.l 20151 ). However, measuring or even identifying 
substructures in lensing clusters has so far been possible 
only under the assumption that light traces mass. Without 
this assumption, the mass function does not seem a viable 
statistic with lensing, especially since lensing gives informa¬ 
tion about the sky-projected mass and projection may wash 
out substructures. 

In this work, we propose a different strategy, which is to 
use a two-dimensional power spectrum as a basis for compar¬ 
ing lensing clusters with simulations. In section [2] we define 
such a power spectrum, which is normalised differently from 
the usual cosmological power spectrum, and has dimensions 
of mass-squared. In section [3] we briefly describe GRALE, 
a free-form lens-reconstruction technique and code. In sec¬ 


tion[4]we apply GRALE to reconstruct the six clusters of the 
Hubble Frontier Field from strong-lensing data, and calcu¬ 
late their power spectra. Then in section [a] we compare the 
clusters with each other and with ACDM simulations, both 
dark-matter only and hydrodynamical. Finally, section[6]has 
general discussion and suggestions for future work. 


2 SUBSTRUCTURE POWER SPECTRUM 

Since every cluster is different, comparison of mass maps 
of observed and simulated clusters is necessarily statistical. 
Properties like number of galaxies, ID density profiles, con¬ 
centration, mass function of the structures/substructures, 
temperature profiles etc. are useful. However, none of these 
quantifies the clustering properties of the halo/cluster, 
which contains important information about its merging his¬ 
tory and evolution. So it is interesting to study statistically 
the cluster mass distribution. The simplest statistic is a two- 
point function. Cosmologically, a two-point function gives an 
excess probability that a local density peak exists close to 
another peak, as a function of the separation between the 
two. It can be well studied in the Fourier space as the power 
spectrum. Generally, a power spectrum P(k) of a 3 dimen¬ 
sional matter density field is defined as: 

P(k) = (\~5(k)\ 2 ) (1) 

where, S(k) is the Fourier transform of the over-density at 
position corresponding to the wave-vector k. 

In this work, we measure the 2D power spectrum from 
the projected mass distribution of the halos: 

Pm(k) = (\M 2D (k)\ 2 ) (2) 

where, M 2 d ( k ) is the Fourier transform of the 2D projected 
mass element at position corresponding to the wave-vector 
k. 

M(k) = J E (x)e^' x d 2 x (3) 

This form also gives a natural normalisation of the function 
to be the square of the total mass of the halo. 

This is not the only way to describeUhe statistica l prop¬ 
erties of a density peak, For example. iHezaveh et al.l ((20141 ). 
using the halo model approach, define 

P2u(k) = J ^\k{k)\ 2 dM , (4) 

where P 2 n(k) is the so-called one-halo term in two dimen- 

dfl 

sions, is the differential mass function (number density 
of halos per unit mass) and k(k) is the Fourier transform of 
the convergence (the normalised surface mass density). The 
framework is based on the assumptions of virialised halos, 
and a functional form for its mass function and radial den¬ 
sity profile of each halo. The total contribution to the power 
spectrum comes by adding the correlation between differ¬ 
ent halos (the two-halo term) and the correlation between 
matter within the same halo (the one-halo term). Within 
a halo, there are also many substructures. If one wants to 
study the clustering properties of matter within a gravita¬ 
tionally bound system like a galaxy cluster (or a merging 
system of galaxies), the two-halo term can be neglected as 
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all the structures are moving within the gravitational poten¬ 
tial of the host and are indifferent to each other’s gravity. 
However, the one-halo term still exists, and at large scales it 
contributes to the Poisson distribution of the substructures 
which drops at the size of the largest substructure. However, 
this can only apply if the cluster (or the halo) is virialised 
so that a smooth density profile can be subtracted in order 
to see the correlation between the residual field. Now, this 
is not a good assumption overall, especially during merg¬ 
ing. As most of the massive galaxy cluster lenses at redshift 
> 0.3 —0.4 are not virialised systems, this form of the power 
spectrum is not intuitive. Therefore, for a general system, 
including ongoing mergers, recent mergers and virialised ha¬ 
los, the correlation between the distribution of matter inside 
the halos must be studied without such assumptions. 

With our definition, if we have a field enclosing only 
one halo with its substructures, and we compute Pm(&) at 
scales larger than the size of the halo, for all those scales one 
gets the same mass and hence, the power above all those 
scales will be a constant, which mimics the Poisson noise 
of the substructures in the halo model picture. Further this 
function drops at scales smaller than the size of the halo, 
which is the largest structure in the field, consistent with 
the halo model intuition. For a virialised halo, this function 
will be smooth and a power law for all scales smaller than 
the size of the halo. However, if the halo is not spherically 
symmetric and in a merging stage the fluctuations in the 
function resemble variation in its power at various scales. If 
one identifies a large ensemble of halos (let’s say in simu¬ 
lations) of nearly the same mass and the same epoch and 
averages all the power spectra, one will get an unbiased trend 
of the clustering properties of such halos statistically. How¬ 
ever, studying individual systems is like studying a random 
realisation of a halo which underlies a mean power spectrum 
and the recovered power spectrum from such a halo should 
be statistically consistent with the mean. 

To ensure this feature, we performed a test. We gener¬ 
ated three density fields which are the random realisation of 
the same power spectrum (oc k~ 3 at small scales) as shown 
in figure |T| The three density grids enclose structures and 
some substructures and look very different morphologically, 
but as they are generated from the statistics, they should 
show a similar power spectrum. In the bottom-right of fig¬ 
ure [l] we show the ratio of the recovered power spectrum 
from each of the grids to the input power spectrum. The 
slope of the recovered power spectrum is calculated by fit¬ 
ting the power law towards the larger values of k, between 
0.68 to 4 (in arbitrary units). All power spectra give nearly 
the same slope as the input one with an error of about 10- 
20%. This test shows that given a density field, the correct 
mean statistic can be obtained within a reasonably small 
error. 

For a pure noise field, we expect same power at all 
scales. We performed this test as well (but we don’t show the 
result in a figure form), and verified that the power spectrum 
is flat. 

The fluctuations in mass distribution at various scales 
can be directly inferred as the presence/absence of substruc¬ 
tures. Larger power at small scale indicates the presence 
of local density peaks, however, a flat power at all scales 
indicates similar structures at all scales or just the noise. 
Also, a smoother distribution of the matter gives rise to 


smaller power whereas sharp peaks correspond to larger 
power. Therefore, we expect the power at small scales to in¬ 
crease as the merger progresses and reach the highest power 
when the system is virialised. In other words, one should 
expect the power on small scales to decrease with redshift. 
This form of the power spectrum can be used in order to 
compare halos statistically, within simulations or observa¬ 
tions or across the two. 


3 LENS RECONSTRUCTION METHOD 

In GRALE llLiesenborgs et aj1l2007l . l2008h the mass distri¬ 
bution is free-form and consists of a superposition of a large 
number of adjustable components (several hundred Plum¬ 
mer lenses). The distribution of these Plummer lenses is 
adaptively determined by GRALE using a multi-objective 
genetic algorithm, to optimise two fitness measures. 

(i) The overlap fitness quantifies the fractional overlap of 
the projected images of the same source back onto the source 
plane. If all images of a source back-project to exactly the 
same area on the source plane, the source fitness for that 
image system is perfect. More generally, the larger overlap 
the better the fitness. It is important to use the fractional 
overlap, as otherwise the fitness measure would be biased 
towards extreme magnifications. 

(ii) The null fitness is a penalty for any spurious images 
implied by a model where none are present in the obser¬ 
vations. A null space is created at the image plane, subdi¬ 
vided into a number of triangles, and the trial solution un¬ 
der study is used to project these triangles onto the source 
plane. Then, the amount of overlap between each triangle 
and the current estimate of the source shape is calculated 
and used to construct a null space fitness measure which is 
the penalty. The penalty is applied only in regions where 
the data make it clear that no images are present; extra 
images are allowed in regions that are difficult to observe 
because, for example, of the presence of nearby bright galax¬ 
ies. Among lens-modelling techniques, GRALE is unique in 
being able to exploit the absence of images as useful infor¬ 
mation. 

Further fitness measures can be defined and in- 
corporated into GRALE, in particular, ti me delays 
llLiesenborgs et ahll2009l : iMohammed et al.ll201ol ). No infor¬ 
mation about the light from the cluster is used — mass- 
traces-light assumptions are completely absent. 

While the genetic algorithm optimizes the properties 
and placement of the component Plummer lenses, the user 
still needs to specify a range for the allowed number of com¬ 
ponents. This is, in effect, an overall resolution of the mass 
map. To set this effective resolution, we run the mass re¬ 
construction process first at coarse resolution, then finer, 
and then coarse again, and choos e an optimal trade-off be - 
tween fitness and resolution — sec Moham med et ahl (l2014h 
for details and tests of this strategy. Finally, the entire cy¬ 
cle is repeated 30 times (while using pre-HFF data) or 40 
times (while using post-HFF data), to generate an ensem¬ 
ble of mass reconstructions. In the remainder of the article, 
we present the average mass map of 30 (or 40) mass re¬ 
constructions for each cluster lens, unless stated otherwise. 
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Figure 1. The first three plots (in the top row: Gridl, Grid2, Grid3) are the three density map realizations generated from the same 
power spectrum with k ~ 3 towards large k. The bottom-left panel shows the actual power spectra of each of the grids along with the 
true power spectrum with slope = -3. The bottom-right panel shows the comparison of the recovered power spectrum from each density 
grid to the input power. These maps were only produced to test the numerical computation of t he power spectrum , and have no central 
cluster. While it is possible to generate random fields constrained to have certain features (cf. iHoffman fe Ribakll 19911 ) these were not 
needed for this simple test. 


The power spectra for each of the reconstructions were com¬ 
puted, and both average power spectra as well as the power 
spectra of the average mass map are presented. For the true 
error estimates, the former is more useful. 


post-HFF data, for clusters MACS0416 and Abell 2744. Ta- 
blc[l]shows the necessary information about the lensing data 
in each of the HFF clusters. In the following figures, the mass 
maps and the corresponding power spectra are marked as vl 
and v3 indicating pre-HFF and post-HFF data respectively. 


4 HUBBLE FRONTIER FIELDS 


4.1 MACS J0416.1-2403 


Hubble Frontier Fields (HFF) survey (PI: J.Lotz, HST 
13498) is a three year Director’s Discretionary Time pro¬ 
gram that devotes a total of 840 orbits to six galaxy clus¬ 
ters plus the accompanying parallel fields. Each field is ob¬ 
served in three HST optical bands (ACS F435W, F606W 
and F814W) and four infra-red bands (WFC3-IR F105W, 
F125W, F140W, and F160W). It is the single most ambi¬ 
tious commitment of HST resources to the exploration of the 
distant Universe through the power of gravitational lensing 
by massive galaxy clusters. All six clusters are early or in¬ 
termediate stage mergers at 2 = 0.3-0.55, with significant 
elongation and hence a non-trivial mass distribution. Each 
had about 10-20 multiply imaged background sources dis¬ 
covered with pre-HFF HST data. Five independent teams 
were tasked with making mass and magnification maps for 
these six clusters. In the analysis of the present paper, we 
use mass maps presented by our team, made with pre-HFF 
lensing data on six clusters as well as two maps made with 


MACS ,10416.1-2403 (MACS0416 henceforth) is a strong 
gravitational lens at redshift of nearly 0.396. It is an elon¬ 
gated galaxy cluster, most proba bly a recent merger or 
a pre-merger (iQgrean et al.l 1201511 , and hence has a non¬ 
trivial mass distribution. It wa s first i dentified by the MAs- 
sive Cluster Survey (MACS; lEbeling et al.l hoiol l. Based 
on its double- peaked surface brightness in the X-rays 
dMann fe Ebelina l2012l l. its recent merger stage is con¬ 
firmed. It is amongst the five high magnification clusters in 
the CLASH (Cl uster L ensing and Su pernovae survey with 
Hubble) project (Postman et al. 20121 1. Detailed mass maps 
were provided bvIZitrin et al.1 d2013ll usin g CLA SH data, and 
Ijauzac et al.l (I2015T1 and ISebesta et al.l (120151 1 using HFF 
data. 

We made two reconstructions of this lens using GRALE, 
one with pre-HFF data (tota l 40 images from 13 sources), us¬ 
ing data from lMerten et aD (1201 il l . as well as data provided 
by J. Richard and D. Coe, and one with HFF data (total 
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Cluster 

Data 

Redshift 

Number of images 

Number of sources 

Redshift range of sources 

MACS-J0416 

pre-HFF 

0.396 

40 

13 

1.82-3.25 

MACS-J0416 

post-HFF 

0.396 

88 

36 

1.00-5.90 

Abell 2744 

pre-HFF 

0.308 

41 

12 

2.00-4.00 

Abell 2744 

post-HFF 

0.308 

55 

18 

2.00-4.00 

Abell 370 

pre-HFF 

0.375 

36 

11 

0.73-3.00 

AS-1063 

pre-HFF 

0.348 

37 

13 

1.22-6.09 

MACS-1149 

pre-HFF 

0.543 

32 

11 

1.23-3.80 

MACS-0717 

pre-HFF 

0.545 

23 

7 

1.80-2.91 


Table 1. Lensing data from different HFF clusters. 
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Figure 2. Left: the density map for MACS0416 cluster using pre-HFF data (top row) and HFF data (bottom row). The crosses mark 
the positions of the lensed images. Right: the highlighted (black) corresponding power spectrum of the cluster on the left with all other 
HFF clusters’ power spectra in the background. 


88 images from 36 sources) where we used images that were 
classified collectively by the HFF teams being of GOLD and 
SILVER quali ty. Many of these images were initially iden¬ 
tified by Ijauzac et al. (120151 ). The two reconstructed mass 
models, vl and v3—with each being an average of 30 and 
40 realisations respectively—are shown in figure [5] (left col¬ 
umn), each with its respective power spectrum (highlighted 
in the right column). Both mass models are very similar at 
large scales, which is also evident from the power spectrum, 


and only differ in the details at small scales. The two mass 
clumps and the elongation can be identified in both models. 
However, the mass models with HFF data (bottom row), 
shows larger power or many detailed substructures at small 
scales. This is expected as HFF data has many more lensed 
images than pre-HFF data and hence provides additional 
lensing constraints at small scales, which leads to increased 
power at larger k’s. 
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A2744-V1 
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Figure 3. Left: the density map for the Abell 2744 cluster using pre-HFF data (top row) and post-HFF data (bottom row). The upper 
left clump could well be an artefact as there are no images there. Right: the highlighted (black) corresponding power spectrum of the 
cluster on the left with all other HFF clusters’ power spectra in the background. 


4.2 Abell 2744 


Abell 2744 is a massive galaxy cluster at redshift 0.308 
and an active merger. It has been studied in various 
wavelengths, for example, it has an exte nded radio halo 


^ i^riovanmm 

(Allen 1998 

et ai. ^ovom et ai. zuuu, A.-ray emission 

Ebeling et ahl 201ol; Kemnner & David 2004; 

Owers et al. 

tified in the 

20111) and various substructures have been iden- 
optical observations (Girardi & Mezzetti 200ll; 

Boschin et £ 

offset betwe 
and baryon 
also a mag 
ground, at i 
ing analyses 

ill 20061). Shan et al. (201()l) show a significant 
en the dark-matter component (from lensing] 
c matter (from X-rays observations). There is 
nified singly imaged supernova in the back- 
edshift 1.35 (Rodney et al. 20151). Various lens- 

have been done (Smail et al. 1997; Allen 

1998 

Cypriano et al. 2004; Merten et al. 2011; Lam et al. 

2014; 


Jauzac et al.ll2014l ). 


We present two mass models of Abell 2744 using 
GRALE, one with pre-HFF data dMerten et al.l l20lJ : 
iRichard et al.l l2014) as well as data provided by J. Richard 
and D. Coe, for a total of 41 images from 12 background 
sources spread over a redshift range of 2 to 4 and another 


with post-HFF data for a total of 55 images from 18 sources. 
Figure [3] shows the reconstructed mass maps and the corre¬ 
sponding power spectra. It shows two very distinct blobs and 
an overall elongated structure, a morphology very likely for 
a major merger which is in fact confirmed by previous stud¬ 
ies. The post-HFF data contains improved identifications 
and positions of the lensing images, and the corresponding 
map (bottom row of figure [3j shows fewer spurious struc¬ 
tures compared to the pre-HFF map (top row of figure [3]|. 

Due to very sharp mass peaks, this cluster, like 
MACS0416, also shows larger power at small scales. 


4.3 Abell 370 

Abell 370 is a strong gravitational lens at redshift 0.375 
and hosts a giant gravitational arc at redshift 0.725. Due 
to its large Einstein radius, it is an ideal target to look for 
high redshift galaxies through its magnification, especially in 
high-density reg ions. There are many published mass mod¬ 
els (for examp le Soucail et al.lll987i : Abdcls alam et al.1ll998l : 
IRichard et al.ll2010l ). We used pre-HFF data from the latter 
work, as well as data provided by D. Coe and A. Koeke- 
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Figure 4. Left: the density map for the Abell 370 cluster using pre-HFF data. Right: the highlighted (black) corresponding power 
spectrum of the cluster on the left with all other HFF clusters’ power spectra in the background. 


moer, to reconstruct its mass distribution. Figure [4] shows 
the resulting mass map. The mass distribution is shallow 
but shows many substructures. The resulting power spec¬ 
trum shows less power on small scales as compared to other 
HFF clusters which is reflected in shallower peaks in its mass 
distribution. 


4.4 Abell S1063 

Abell S1063 is a galaxy clus ter at redshift 0 .348 . We used 
pre-HFF data identified by iRichard et al.l (120141 ); our re¬ 
construction used a total of 37 images from 13 background 
sources. 

Figure [5] shows the reconstructed mass map along 
with the corresponding power spectrum. It shows two mass 
clumps, which lead to larger power at intermediate scales, 
but due to the shallowness of the peaks, it loses power at 
small scales. 


4.5 MACS J1149.5+2223 

MACS J1149+2223 (MACS 1149 hereafter) is an X-ray 


bright cluster at redshift 0.543. It has been 
various authors for its rich strong lensing 1 

studied by 
Smith et al. 

200S|; Zitrin &; Broadhurst 2009 

: Zitrin et al. 

2011], 120151: 

Rau et al. 20lj Johnson et al. 

2014k Sharon & Johnson 


20151). The re is also a multiply imaged supernova 
llKellv et al.l 120141 . 120151) hosted by a face-on spiral galaxy 
at redshift 1.49. 


Figure [6] shows the reconstructed mass model of 
MACS 1149 using GRALE with pre-HFF da ta (total 32 
image s from 11 background sou rces) from ISmith et al.l 
(12009) : IZitrin fc Broadhurstl J2009I ). and data provided by 
M. Limousin. The mass model favours two dominant peaks 
and one sub-peak, elongation and many substructures. The 
peaks are shallower than those in Abell 2744, also expected 
as it has nearly the same mass as Abell 2744 but higher 
redshift. The corresponding power spectrum is shown in the 
right column of the same figure. 


4.6 MACS J0717.5+3745 


MACS J0717.5+3745 (MACS0717 hereafter) is a strong 
gravitational lens at redshift 0.55, classifie d as the most dra¬ 
matic merger in X-ray/optical analysis bv lMann fe Ebelind 
( 120121 ) . 

We used pre-HF F da ta from iLimousin et ah (120121) : 
IRichard et al.l d2014l) and iMedezinski et ah ( 2015 ) to re¬ 
construct its mass distribution using GRALE. Figure [7] 
shows the resulting mass map and the corresponding power 
spectrum. Due to the clear mass peaks and substructures, 
the power spectrum shows increased power at intermediate 
scales, however, due to very shallow peaks the power spec¬ 
trum at small scales drops and is amongst the lowest of the 
HFF clusters. 


5 COMPARING THE CLUSTERS 
5.1 Comparing HFF clusters 

We can now compare the clustering properties of HFF clus¬ 
ters using the lensing power spectrum Pjyi(fc) • The recon¬ 
structed mass distributions of the HFF clusters are mor¬ 
phologically different but show similar statistical structures 
which is very much evident from comparing the respective 
power spectra. Figure [H] shows all the power spectra for 
HFF clusters in thick lines along with their sample vari¬ 
ance, which represents scatter across all available modes for 
a given k. 

Let us first consider two reconstructions for MACS0416, 
using pre-HFF and HFF data. Both mass maps are very 
similar at large scales but differ in detail at small scales. 
The same effect can be noticed in the power spectrum (see 
figure[8]): up to k ~ 3 x 10 _5 pc~ 1 both look nearly the same, 
but at smaller scales the mass map reconstructed with pre- 
HFF data starts to lose power. This is a consequence of 
the fact that HFF data has over two times as many lensed 
images as pre-HFF data (see table |TJ, allowing for higher 
resolution in the reconstructed mass maps. This leads to a 
higher power at small scales, where the difference in the two 
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AS1063-V1 




Figure 5. Left: the density map for the AS1063 cluster using pre-HFF data. Right: the highlighted (black) corresponding power spectrum 
of the cluster on the left with all other HFF clusters’ power spectra in the background. Note that there are no images in the upper right 
corner of the map, so the massive clump in that region could well be an artefact. 




Figure 6. Left: the density map for the MACS1149 cluster using pre-HFF data. Right: the highlighted (black) corresponding power 
spectrum of the cluster on the left with all other HFF clusters’ power spectra in the background. There are no images beyond X=25”, 
so the two clumps on the right could well be artefacts. 


power spectra is almost an order of magnitude. A similar 
trend can be seen in the two versions of Abell 2744. 

Power spectra for Abell 370 and MACS 1149 look very 
similar to each other. If we look at the corresponding 
mass maps, both show two dominant peaks with similar 
width, shallow as compared to the surroundings with similar 
smoothing. Similar arguments can be made for Abell 2744 
and MACS0416-II clusters. The clusters show similar power 
at small scales which is expected from the very steep mass 
peaks in the two clusters. 

It is important to note that the smallest scale (largest 
k) to be trusted depends on the density of images. For 
example, MACS0416-v3 has 4 times more images than 
MACS0717, (assuming roughly the same area for both in 
terms of arcsec 2 ), so the typical linear image spacing is 2 
times larger in MACS0717, leading to poorer mass resolu¬ 
tion. And, in fact, looking at figure [8j the k value where 


Pm(^) of MACS0717 begins to drop (k ~ 2 x 10 5 pc x ) is 
about a factor of ~ 2 smaller compared to that of MACS0416 
(k ~ 4 x 10 -5 pc -1 ). 

There is also a tendency for low redshift clusters to have 
higher power at small scales. Thus, Abell 2744 reconstructed 
with 41 images has larger power than MACS0416-vl (40 im¬ 
ages) at small scales. Similarly, MACS 1149 has less power at 
small scales than A370 and AS1063, again with similar num¬ 
bers of lensed images. This is also intuitive as low redshift 
clusters have had more time to build small-scale substruc¬ 
tures. Therefore, the power at small scales may be indicative 
of the age of the cluster. 

5.2 ACDM simulation clusters 

We used 22 clusters of galaxies from dark-matter only 
(DMO) simulations and an additional 22 from hydrodynam- 
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MACS0717-V1 


10 



X [arcseconds] 



Figure 7. Left: the density map for the MACS0717 cluster using pre-HFF data. The upper left and bottom right clumps could well be 
artefacts as there are no images there. Right: the highlighted (black) corresponding power spectrum of the cluster on the left with all 
other HFF clusters’ power spectra in the background. 




DM0 
AGN 
A2744—vl 
A370 — vl 
AS1063 —vl 
MACS0416 — vl 
MACS0717 — vl 
MACS 1149 — vl 
MACS0416 - v3 
A2744 —v3 


DM0 
AGN 
A2744—vl 
A370 — vl 
AS1063 —vl 
MACS0416 — vl 
MACS0717 — vl 
MACS 1149 — vl 
MACS0416 - v3 
A2744 —v3 


Figure 8. All power spectra: clusters in thick solid lines and simulations in thin solid lines. Left: showing power spectra calculated for 
the average mass map of each cluster. The error bars are the scatter in different modes at a given k (or the sample variance). Right: 
each power spectrum is the average of 30 (40 in case of v3) different power spectra calculated for each mass map in the ensemble for the 
respective cluster. The error bars are the standard deviation of the power spectra, including contributions from sample variance, as well 
as the scatter between different reconstructions of the same cluster. 


ical simulations which i nclude AGN feedback , in the mass 
rang e l-3x 10* 4 Mq from lMartizzi et alJ (1 20141 ) using RAM¬ 
SES (|Tevssiei]|2002 ). The mass range is chosen to be compa¬ 
rable to that of the HFF clusters but the simulated clusters 
are all at redshift zero. We then projected them along the 
three axes, for a total of 66 projected clusters in each case. 
Projection does introduce some correlation in the sample, 
but we assume such correlations in the sample, but we ne¬ 
glect this effect. In contrast to HFF clusters, the simulation 
clusters are virialised but a few show more than one core. 


The measured power spectra are shown as thin lines in figure 

m 


We also measured the mean power spectrum in each 
case, DMO and AGN, which is shown in black in the left 
panel of figure^ The error bars show the standard deviation 
of Pm(&), including contributions from scatter in different 
modes at a given k as well as the scatter between different 
clusters. In the right panel of figure 0 we show the rela¬ 
tive difference between the mean power in DMO and AGN 
simulated clusters. Comparing the mean power spectra of 
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Figure 9. Left: comparing the mean power spectrum for dark-matter only simulations (DMO), baryonic simulation (AGN) and different 
reconstructions using GRALE and LENSTOOL for MACS0416-v3 (using HFF data only). The vertical lines show the 10 and 1 arc second 
scales at the cluster’s redshift z = 0.396. Right: the relative difference between the mean power spectra of DMO and AGN simulations. 



k[pc *] 


Figure 10. Comparing average power spectrum of HFF clusters to that of the simulated clusters. 


DMO and AGN, the latter shows a deficit in power at large 
scales and a boost at smallest scales. The deficit is nearly 
10%. The AGN clusters have the influence of an AGN at the 
centre that drives the gas outside the cluster with its feed¬ 
back process. Losing this mass results in the deficit of power 
at those scales. However, these AGN clusters contain a cen¬ 


tral stellar component in the form of the brightest cluster 
galaxy (BCG), which increases the mass at the very centre 
of the cluster resulting in a boost in the power at smallest 
scales. The trend is systematically consistent with previous 
findings. 
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5.3 Comparison between HFF clusters and 
ACDM simulation clusters 

In figure[8] we show Pm(&) for all six HFF clusters as well as 
both DMO and AGN for all the simulated clusters. The lens- 
ing clusters systematically show less power at small scales 
compared to the simulated ones. We discuss the possible 
reasons for this in section [6] 

5.4 Comparison with LTM models 

In the left panel of figure 0 we compare the power spectrum 
measured from the mass maps reconstructed using GRALE 
to that made using other methods: Sharon et al and the 
CATS groupQ] Bo th use the LTM method LENSTOOL of 
iJullo et al.l (|2007T ) as the reconstruction technique. Both of 
these reconstructions show much larger power at small scales 
compared to the GRALE reconstruction. This is the result 
of the fact that LENSTOOL uses information from the in¬ 
dividual cluster galaxies of the lens, and hence has much 
steeper density gradients at small scales which results in 
larger power. This effect is also seen in the galaxy-mass cor¬ 
relation function of M ACS04 16 which shows a sharp peak at 
zero separation (ISebesta et al.ll2015i h The power spectrum 
of parametrically reconstructed clusters is a mixture of pri¬ 
ors and lensing information, and when the number of images 
is small, priors are the dominant contribution. In contrast, 
free-form methods like GRALE base the mass distribution 
on the lensed images alone, without relying on visible galax¬ 
ies. 


6 CONCLUSIONS AND DISCUSSIONS 

In this article, we present free-form mass models of six Hub¬ 
ble Frontier Field clusters using strong gravitational lens¬ 
ing pre-HFF data. We also studied two of the six clusters, 
MACS0416 and Abell 2744, using both pre-HFF and HFF 
data (with many more lensed images). For each cluster, a 
total of 30 mass maps (40 for v3 versions) are made using 
a non-parametric lens inversion technique called GRALE, 
and the average mass map is presented. All the mass maps 
show elongation, multiple cores and many substructures in 
the mass distribution implying a recent major merger. The 
lensing data from recent HST-HFF observations is rich in 
lensing images allowing a very precise identification of sub¬ 
structures to be done with greater confidence. These mass 
maps can be used to study various aspects of the structure 
formation and merging stages of large collapsed clusters of 
galaxies. 

We measured the power spectra of these sky-projected 
mass maps (30 for each cluster), and presented both av¬ 
erage power spectra, as well as the power spectra of the 
average mass map for each cluster (figure [8|. The error es¬ 
timates of the power spectra are either given as the scatter 
in different modes at a given k, also known as sample vari¬ 
ance, or a combination of the sample variance and the stan¬ 
dard deviation of the power spectra for the same cluster. 
Power spectra give a statistical description of the clustering 
properties of the mass distribution and encode information 

1 https://archive.stsci.edu/prepds/frontier/lensmodels/ 


about the abundance of substructures and their contrast 
with the background. There are tentative indications that 
low redshift clusters have systematically higher power at 
small scales as compared to high redshift clusters, presum¬ 
ably because low redshift clusters are closer to their virial 
equilibrium than high redshift clusters. We argue—and il¬ 
lustrate using pre-HFF and HFF maps of MACS0416—that 
using a larger number of lensed images in a cluster leads to 
a better-constrained mass distribution, and that the power 
spectrum can be recovered up to much smaller scales putting 
stronger constraints on our understanding of the substruc¬ 
tures. On the other hand, parametric methods that explic¬ 
itly include cluster galaxies in the models have much larger 
power on small scales even with fewer lensed images indi¬ 
cating that the mass maps are dominated by the priors, 
especially at small scales. 

We compared the power spectra of the HFF clusters 
with those from ACDM dark-matter only and hydrodynam- 
ical simulations at redshift zero. Figure [10] summarizes our 
comparison between HFF clusters and ACDM simulation 
clusters. We average the power spectra of the six HFF clus¬ 
ters and compare them with the average power spectra of 
the simulated clusters. The average power of the clusters is 
systematically steeper than that of the simulated clusters 
and hence have less power at small scales. This may be due 
to one or more of the following reasons: 

• The redshift of the HFF clusters are higher than the 
simulated clusters. We are comparing non-virialised HFF 
clusters with the virialised halos of the simulations. At later 
times, the mass distribution becomes more and more clumpy 
and the substructures pull more mass from the background, 
which results in a higher power at small scales. Therefore, it 
is possible that the lensing clusters in the local Universe may 
show similar power as the simulated clusters at redshift zero. 
Conversely, this can be checked by comparing the power 
spectrum of HFF clusters with that of simulated clusters of 
similar masses in the redshift range 0.3-0.5. 

Another effect of the redshift difference is that the angular 
resolution is lower for the higher redshift clusters. However, 
this effect is roughly 20%, too small to account for the ob¬ 
served differences in the power spectrum. 

• A second possible reason relates to the data and method 
we are using to reconstruct the clusters. A larger number 
of lensed images gives additional constraints on the mass 
distribution, and this increased resolution leads to a boost 
in the power spectrum. This can be seen in figure [8] where 
we show the power spectrum of MACS0416 using both pre- 
HFF and HFF data: HFF data shows larger power at small 
scales as compared to pre-HFF data. 

Figure [2] shows that the power at small scales also de¬ 
pends on the assumption of the reconstruction method. It is 
possible that GRALE does not have enough resolution and 
hence loses power at smallest scales. On the other hand, 
LENSTOOL maps seem to have much more power than ex¬ 
pected at those scales and are possibly dominated by priors. 

• Finally there is the possibility that the simulations are 
lacking some physics that needs to be taken into account in 
order to simulate realistic halos. 

The second item above can be tested through a more ex¬ 
tended pipeline, in which clusters are generated from state- 
of-the-art simulations, lensing data are generated from them, 
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and then analysed by independent groups using different 
techniques. Such tests are currently in progress. It will be in¬ 
teresting to measure the power spectrum of the original mass 
distribution and the reconstructed one using different lens 
inversion methods. The current analysis expects that the 
power spectrum based on GRALE maps will match the true 
power spectrum at large scales and ultimately lose power 
depending on the density of lensed images, whereas LTM- 
based methods may continue having larger power even when 
the true power drops. We leave this analysis for future work. 

Independent of lensing, it will also be interesting to 
study the evolution of the power spectrum as a probe of sub¬ 
structures and the merging history of large collapsed objects 
in simulations. For example, a cosmological simulation can 
be set up with a volume large enough to produce 20-50 clus¬ 
ters of galaxies, in a narrow mass range. The power spectra 
can be calculated for all the clusters in the mass range at 
different redshifts, and the evolution of the average power 
spectrum can be studied. This will give us insight into the 
evolution of clustering properties in a merging system. We 
leave this analysis to future work as well. 
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